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Abstract 

Organic soils are an important source of N2O, but global estimates of these fluxes remain uncertain because measurements 
are sparse. We tested the hypothesis that N2O fluxes can be predicted from estimates of mineral nitrogen input, calculated 
from readily-available measurements of CO2 flux and soil C/N ratio. From studies of organic soils throughout the world, we 
compiled a data set of annual CO2 and N2O fluxes which were measured concurrently. The input of soil mineral nitrogen in 
these studies was estimated from applied fertilizer nitrogen and organic nitrogen mineralization. The latter was calculated 
by dividing the rate of soil heterotrophic respiration by soil C/N ratio. This index of mineral nitrogen input explained up to 
69% of the overall variability of N2O fluxes, whereas CO2 flux or soil C/N ratio alone explained only 49% and 36% of the 
variability, respectively. Including water table level in the model, along with mineral nitrogen input, further improved the 
model with the explanatory proportion of variability in N2O flux increasing to 75%. Unlike grassland or cropland soils, forest 
soils were evidently nitrogen-limited, so water table level had no significant effect on N2O flux. Our proposed approach, 
which uses the product of soil-derived CO2 flux and the inverse of soil C/N ratio as a proxy for nitrogen mineralization, 
shows promise for estimating regional or global N2O fluxes from organic soils, although some further enhancements may 
be warranted. 
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Introduction 

Although organic soils occupy only 3 % of the Earth's land area, 
they contain approximately 40% (610 Pg) of the terrestrial soil 
organic carbon (SOC) [1]. Climate warming and human 
disturbance such as drainage and cultivation are expected to 
accelerate carbon decomposition in organic soils, and the 
decomposition of SOC can facilitate the release of mineral 
nitrogen which can then be utilized by denitrifying and nitrifying 
bacteria to produce the potent greenhouse gas NgO [2,3]. N2O 
emissions from organic soils under agricultural use in Nordic 
countries were on average four times higher than those from 
mineral soils, indicating that N2O derived from SOC decompo- 
sition dominates overall fluxes [4]. However, no consistent and 
quantitative relationship has been reported for N2O emission and 
organic carbon decomposition in organic soils. 

Organic carbon and nitrogen in soils, plant and microbial 
biomass are usually covalently bonded at relatively constant ratios. 
It is thus logical to expect that N2O and CO2 originated from 
SOC decomposition should be closely linked. Some studies have 
indeed found a significant relationship between soil N2O and CO2 
emissions at the site level [5,6]. This relationship, however, was 
weaker when data were pooled across sites or ecosystems [7, 8]. The 



variability of soil C/N ratio may be one of the important factors 
undermining the correlation for organic soils. The C/N ratio in 
organic soils ranges from 50~100 in weakly decomposed peat to 
12~35 in highly decomposed peat [9]. The supply of mineral 
nitrogen from SOC decomposition is the outcome of two 
concurrent and oppositely directed microbial processes - nitrogen 
mineralization and immobilization [10]. Soils with a high C/N 
ratio may be characterized by rapid immobilization of nitrogen 
and soils with a low C/N ratio by higher net nitrogen 
mineralization and a surplus of available NH^"*" and NO3 [1 1]. 
A negative relationship has accordingly been shown for C/N ratio 
of soils and N2O fluxes [9]. Similar to the relationship between 
N2O and CO2 emissions, the correlation of N2O emission with soil 
C/N ratio tended to be weak when the data from dilferent sites at 
larger scales were included [4, 1 2] , which makes it difficult to scale 
up N2O fluxes by CO2 emissions or C/N ratio alone from 
individual sites to regional scales. In view of the coupling of soil 
carbon and nitrogen processes and the bridging function of C/N 
ratio, we hypothesized that a combination of soil CO2 emission 
and C/N ratio would likely provide better measurements of N2O 
emission at larger scales. In fact, Mu et al. [13] have linked N2O 
flux to soil mineral nitrogen as estimated by CO2 emission and C/ 
N ratio for agricultural mineral soils. To our knowledge, no such 
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kind of attempt has ever been made for organic soils. The aim of 
this study was therefore to determine: 1) if NgO flux from organic 
soils is related to soil mineral nitrogen input estimated from 
heterotrophic respiration divided by soil C/N ratio (a derived 
measure of soil nitrogen mineralization) plus fertilizer nitrogen; 
and 2) whether or not the relationship is sufficiently robust to serve 
as an approach for estimating NgO flux from organic soils. 

Materials and Methods 

Data source 

To test the hypothesis, we collected journal-published data of 
N2O and CO2 emissions measured simultaneously in the fields on 
peadands or histosols for which the carbon and nitrogen content 
or ratio of the organic matter in the upper layers of the soil has 
been reported. Occasional and short-period flux measurements 
were not used and only data on annual emissions were considered. 
For long-term measurements, we used annual estimates rather 
than multi-year averages to reflect temporal variability. Annual 
emissions were directly reported by authors or estimated from 
points in the figures of publications. The final dataset comprised of 
122 field measurements from 28 geographical sites (Table SI). Of 
all data, only 12 measurements at 9 sites were from the tropical 
regions and the rest were from the temperate regions. Most of the 
flux measurements were made using closed chamber technique 
with samphng frequency varying from 1-3 times per week to once 
per month. Other factors such as soil pH and water table level, if 
reported, were also recorded in the database. Readers should refer 
to the original papers for a more complete presentation of the 
data. 

Estimation of soil mineral nitrogen input 

The CO2 emission measured in bare soils can be taken as the 
proxy of SOC decomposition or heterotrophic respiration [14]. 
There are limited studies in which CO2 emission was measured in 
bare soils (Table SI). For the COj emissions measured in soils with 
plants, the contribution of heterotrophic respiration or SOC 
decomposition was estimated using the following equation adapted 
from Bond-Lamberty and Thomson [15]: 

Rh = 10eP-22 + 0-87M«(/10)] (1) 

where Rh is heterotrophic respiration and R, is total soil 
respiration (kg C ha ' yr '). 

The nitrogen mineralization rate from soil organic matter was 
then calculated using the following equation [13]: 

Nm=RhlScN (2) 

where Nm is the gross nitrogen mineralization (kg N ha~' yr~') 
and ScN is soil C/N ratio. 

The mineralized nitrogen from soil organic matter decompo- 
sition and the inorganic nitrogen from chemical fertilizers 
constitute the total input of soil mineral nitrogen (N^f). Atmo- 
spheric nitrogen deposition, as another important external source 
of soil mineral nitrogen, was not considered for our study since 
there were few papers reporting it. 

Statistical analysis 

The dataset in the current study is of unbalanced nature with 
observations collected from peer-reviewed papers rather than from 
systematically designed experiments. Accordingly, the effects of 
soil mineral nitrogen input and other variables on NjO flux were 



analyzed using the mixed model-REML estimation method of 
SAS/MIXED procedure (version 9.3), which is suitable for 
handling unbalanced data. The values of NjO flux were first 
natural-log transformed to normalize their distribution and then 
analyzed by the following model: 

ln(fN2o) = constant -I- ln(Nmf ) -I- pH -I- WT -I- NS,- -I- E cosysj 
-I- NSi X ln(Nmf ) + E cosysj x ln(Nmf ) + E cos ysj x WT 

where fN2o is the NjO flux; Nn^f, pH, WT, NSi and Ecosysj are the 
fixed effects of mineral nitrogen input, soil pH, water table level, 
nitrogen source (i is mineralized nitrogen only or a combination of 
mineralized nitrogen and inorganic nitrogen from chemical 
fertilizers), and ecosystem type (j is forest or non-forest type), 
respectively. A preliminary check of the data showed that the 
general trend of N2O flux in forest system differed from grass and 
cropland, so the ecosystems were simply classified into two 
subclasses as forest and non-forest. Some two-factor interactions 
were also included in the model. A significant level of ^ = 0.05 was 
used to determine if a given variable or interactive effect was kept 
in the model to further seek solutions for fixed effects. Four 
negative values of NjO flux reported by Inubushi et al. [16] and 
Mojeremane et al. [17] can not be subjected to log-transformation 
and ^\•(■rl' not included in the analysis. In addition to determination 
coefficient (i.e., R^ value), concordance between observed N2O 
fluxes and model fits was also analyzed using Lin's concordance 
correlation coefficient (CCC, Stata SE 12.0) to assess the goodness- 
of-fit of the finalized models. The resulting CCC was interpreted 
using the benchmarks described by Klevens et al. [18] as follows: 
<0.20 is considered virtually no agreement; 0.21-0.40 is 
considered slight; 0.41-0.60 is considered fair; 0.61—0.80 is 
considered moderate; and 0.81-0.99 is substantial. 

Results 

As shown in Table 1, soil pH, soil mineral nitrogen source (NS) 
and ecosystem type did not affect the annual N2O flux (/;>(). 05), 
while the input of soil mineral nitrogen (Nj„|) and water table level 
(WT) had significant effects on N2O flux {p<Q.Q\). The F value of 
Nmf was the biggest, indicating the input of soil mineral nitrogen 
was the main factor controlling N2O emission in organic soils. The 
two-factor interactive effects between NS, N,-jjf, WT and ecosystem 
type on N2O flux were not statistically significant (/)>0.05). 

Only the significant variables were then kept in the model to 
solve the estimates for their effects. Two models with different 
combinations of independent variables are shown in Table 2. The 
first model was the simplest one with N^f as the single independent 
variable. The second model was expanded by adding the effect of 
water table level. The 95% confidence intervals of the estimated 
effect of Njnf were overlapped for different models. The models 
indicated that N2O flux was positively correlated with N^r and 
negatively with water table level. Using the estimated effects and 
the variables in the dataset allowed a comparison between 
predicted and observed annual N2O fluxes from organic soils. 
The variable N^^f explained up to 69% of the variability in the 
overall data of observed N2O fluxes (Fig. 1), while the addition of 
water table level increased the explanatory ability to 75% (Fig. 2). 
When the overall data were further divided by ecosystem types, 
the performance of models was somewhat different (Fig. 1 & 2). 
For forest, the determination coefficient (R^) was nearly stable at 
the value of 0.63 for both models. In contrast, the introduction of 
water table level into models slightiy improved the fitted results for 



PLOS ONE I www.plosone.org 



2 



May 2014 I Volume 9 | Issue 5 | e96572 



Linking N'^O Emission to Soil IVIineral N Input 



Table 1. Results of type III tests of fixed effects. 





Effect 


Numerator DF 


Denominator DF 


F Value 


Pr>F 



Nmf 1 


96 


13.16 


0.0005 


pH 1 


96 


1.43 


0.2344 


m 1 


96 


5.15 


0.0255 


NS 1 


96 


0.10 


0.7472 


Ecosystem 1 


96 


0.70 


0.4040 


NSxN„, 1 


96 


0.11 


0.7426 


WTxN,„f 1 


96 


3.20 


0.0767 


Ecosystem xNmf 1 


96 


0.21 


0.6506 


Ecosystem xWT 1 


96 


2.17 


0.1437 



Nmf, the mineral nitrogen input to soil; WT, water table level; NS, the source of soil mineral nitrogen. 
doi:l 0.1 371 /journal.pone.0096572.t001 



non-forest systems with values increasing from 0.59 to 0.69. 
This indicated that tlie input of mineral nitrogen was the most 
important predictor of N2O flux, while water table level was a 
weak predictor of N2O flux and appeared to be dependent on 
ecosystem type. 

The slope of regression lines in Fig. 1 & 2 ranged from 0.50 to 
0.75, indicating that the relationship strays from the ideal 1:1 line. 
Therefore the concordance correlation coefficient (CCC) between 
observed and predicted N2O fluxes was calculated to measure 
robustness of the models. For the overall data with log- 
transformation, the concordance was substantial with the CCC 
ranging from 0.82 to 0.86 for the two models. When the log- 
transformed data were converted to actual NjO fluxes, however, 
the cluster of fluxes greater than 15.0 kg N ha ' yr ' was found to 
be distincdy underestimated. The CCC for this cluster of data 



ranged from —0.002 to 0.16 and showed virtually no agreement, 
suggesting that some important factors responsible for these high 
fluxes were not accounted for by the models. For the rest of the 
data (103 fluxes out of 1 18), the CCC (ranging from 0.63 to 0.68) 
stiU showed a moderate concordance. 

The variable N^f in the models can be decomposed into soil 
heterotrophic respiration (Rh), C/N ratio and inorganic nitrogen 
rate from chemical fertilizer (Nj). The mixed procedure analysis 
indicated that each of these components of N^^f had a significant 
influence on NjO flux (/)<0.001), with R^ and Nf being positively 
related to N2O flux and C/N ratio negatively related to N2O flux. 
Soil carbon and nitrogen contents, which could replace the 
variable of C/N ratio, were also significantly negatively or 
positively correlated with N2O flux (/;<0.001). The fitting 
efficiency between observed and predicted N2O fluxes by models 
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Figure 1. Correlation of observed fluxes of N2O from organic 
soils and predicted values (kg N2O-N ha^) by model 1 as 
presented in Table 2: IniN^O flux)= 1.8685 ln(Nn,f)-9.0314. 

Solid line shows linear regression fit for the overall data: y = 0.69x+0.13, 
= 0.69. Long-dashed line shows linear regression fit for the non-forest 
system: y = 0.50x4-0.68, R^ = 0.59. Short-dashed line shows linear 
regression fit for the forest system: y = 0.69x-0.33, R'^ = 0.63. 
doi:1 0.1 371 /journal.pone.0096572.g001 
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Figure 2. Correlation of observed fluxes of N2O from organic 
soils and predicted values (kg N2O-N ha^) by model 2 as 
presented in Table 2: ln(N20 flux) = 1.5374 ln(N„,f)-0.0221 
WT— 8.2334. Solid line shows linear regression fit for the 
overall data: y = 0.75x4-0.08, R^ = 0.75. Long-dashed line shows 
linear regression fit for the non-forest system: y = 0.65x4-0.48, R'^ = 0.69. 
Short-dashed line shows linear regression fit for the forest system: 
y = 0.65x-0.35, R^ = 0.63. 
doi:10.1371/journal.pone.0096572.g002 
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using the above-mentioned components of Nmf as inputs were 
nearly the same as those of models using Ni-^f itself (data not 
shown). 

Discussion 

Previous studies have linked NjO flux directly to either CO2 
flux or soil C/N ratio [5,8,9]. In this study, soil CO2 emission and 
C/N ratio were combined to estimate mineral nitrogen input, and 
the latter accounted for up to 69% of the variability of N2O fluxes 
from organic soils with various properties, land management 
practices and climates. Soil CO2 flux or C/N ratio alone explained 
only 49% and 36% of the overall variability of N2O fluxes, 
respectively (Fig. 3). This suggests the necessity of combining soil 
CO2 flux and C/N ratio for predicting N2O flux on a large scale. 
Of course, soil CO2 flux and C/N ratio can be independendy 
incorporated into the same models, but the interpretation of such 
models would be relatively complicated and evasive since there are 
various mechanisms which may explain the control of CO2 flux 
and C/N ratio over N2O flux [8,9,19]. In contrast, the quotient of 
soil CO2 flux and C/N ratio can well represent in theory the gross 
nitrogen mineralization [20] , and the implication of models using 
such a quotient as input is straightforward and self-evident in the 
importance of mineral nitrogen input for regulating soil N2O flux. 
There is no significant diflFerence in the influence of different 
sources of mineral nitrogen on N2O flux (Table 1), suggesting that 
the simplified models might also be suitable for evaluating the 
effect of mineral nitrogen from other sources such as atmospheric 
deposition, though this idea needs further verification. 

A negative relationship between N2O flux and groundwater 
level has been observed for individual sites [21,22], and still holds 
at a large scale as shown in this study. This is logical simply 
because high moisture with increasing water table level can limit 
N2O emission from soils due to the low availability of nitrate and/ 
or efficient reduction of N2O to N2 through denitrification [16,23], 
while the lowering of water table increases oxygen penetration into 
the peat and enhances the decomposition of organic matter, as 
indicated by the negative relationship between heterotrophic 
respiration and water table level (R^ = 0.31, /)<0.0001). It has 
been reported that the control of soil water content or water table 
level over N2O flux is important only when soil is not nitrogen 
limiting [24,25]. In this study, the percentage of observations with 
Nnif greater than 150 kg N ha~' was only 19% for forest, but up to 
87% for non-forest systems (Table SI). This suggests that forest soil 
is nitrogen limiting when compared with non-forest systems, which 
may be responsible for the insensitivity of N2O flux to water table 
level for forest systems (Fig. 1& 2). Besides the input of mineral 
nitrogen, forest differs from non-forest systems in many other 
factors, such as vegetation, below-/above-ground biomass, litter 
fall, soil compaction, and land management practices, all of which 
can influence N2O flux but are not considered here due to limited 
and unsystematic information in literature sources of the current 
dataset. To fill the gap, ecosystem type was used as a proxy 
variable that we tried to incorporate into models; however, 
statistical analysis showed that its effect was not significant 
(Table 1). 

It should be acknowledged that the models described here were 
dependent on simplifying assumptions that can introduce error. 
That is, the gross nitrogen mineralization was estimated from 
carbon mineralization and soil C/N ratio by assuming that the 
rate of carbon mineralization is the same as the rate of respiration 
and the C/N ratio of mineralized organic matter is the same as 
that of the bulk soil organic matter. In fact, carbon and nitrogen 
mineralization from soils originates from decomposable fractions 
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Figure 3. Correlation of observed fluxes of N2O with estimated soil heterotrophic respiration (Rh, left panel) and C/N ratio (right 
panel). 
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of organic matter with different C/N ratios [26]. Most likely, the 
ratio of carbon evolved/ nitrogen mineralized is much wider than 
the bulk soil carbon to nitrogen ratio [27,28]. This indicates that 
gross nitrogen mineralization might be over- or under-estimated if 
bulk soil C/N ratio was used in equation 2. The respiration 
process is also not exactly identical to carbon mineralization. The 
amount of carbon that is ultimately lost through respiration 
depends on how effectively the decomposer community converts 
mineralized carbon to biomass [29]. Similarly, the amount of 
nitrogen that is ultimately available to denitrifier or nitrifier for 
producing N2O depends on how effectively the decomposer 
community converts mineralized nitrogen to biomass and plants 
compete with microbes for mineral nitrogen [10,20]. Empirical 
relationships have been established between nitrogen and carbon 
mineralization in studies performed usually under laboratory 
conditions [30]. Different organic matter fractions or their C/N 
ratios, and varying microbial use efSciency of carbon and nitrogen 
have also been proposed to predict nitrogen release [20,29]. 
However, these relationships are strongly dependent on the 
experimental conditions in which they have been established. 
Moreover, the current dataset is based on the in situ measure- 
ments in the field environment and contains only the basic 
information of respiratory carbon and bulk soil C/N ratio, thus 
necessitating the above-mentioned assumptions to estimate min- 
eralized nitrogen. Such simplifications and assumptions may bring 
uncertainties, but it is necessary in some cases to understand the 
general trends and probabilistic nature of the environment [31]. 

N2O emission from soils is of small magnitude and highly 
variable in space and time, and is thus very difficult to estimate. 
The measurement of soil N2O flux also requires intricate 
techniques along with a lot of time and labor. In contrast, soil 
CO2 emission is controlled primarily by soil temperature and 
moisture, and is relatively easy to measure or predict [32,33]. In 
addition, the estimates of soil respiration are currently more widely 
available than those of soil N2O emission. The models developed 
in this study showed a promising approach to estimating N2O 
emission from organic soils by using soil C/N ratio and CO2 
emission data derived from measurements or biogeochemical 
modeling. It should be mentioned, however, that several aspects of 
the information in the current dataset might impose uncertainties 
on thes(- models. First, soil heterotrophic respiration was simply 
estimated from total soil respiration using a universal relationship 



between them [15], but the relative contribution of organic matter 
decomposition or heterotrophic respiration would vary over time 
and depend on root respiration of the growing plants [8] . Second, 
the majority of the global organic soils are distributed in the boreal 
and sub-arctic regions and about 10%-15% in the tropical 
countries [1,3], but most of the current data came from northern 
Europe, indicating that the models developed in the present study 
might be biased to the temperate regions. 

Conclusion 

A fairly large number of data were collected to explore the 
relationship between annual N2O emission and multiple variables 
for organic soil by a mixed-model analysis, and the input of soil 
mineral nitrogen was found to be the most useful pr<;di( tor for 
N2O flux. Soil mineral nitrogen was supposed to be composed of 
organic nitrogen mineralization as estimated by CO2 emission and 
soil C/N ratio, thus providing a possibility for upscaling N2O 
emission from organic soils by use of regional soil databases 
including information on C/N ratio and carbon storage change or 
CO2 emission data. The approach proposed here may have 
validity as a whole, but needs further evaluation and advancement 
before practical application due to uncertainties associated with 
simplifying assumptions and a regionally unbalanced data source. 
A better understanding of the processes of carbon and nitrogen 
mineralization and their stoichiometric relationship as well as 
additional experimental data from organic soils outside of 
temperate Europe regions wiU help to improve the relationship 
established in this study. 

Supporting information 

Table SI Annual emissions of NjO and CO2 from 
organic soils and estimates of soil mineral nitrogen 
input. 
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